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ABSTRACT 

Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accom¬ 
panied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash 
has been detected. Recent work suggests that the GeV-l-optical flash is emitted by the 
external blast wave from the explosion in a medium loaded with copious pairs. The 
full light curve of the flash is predicted by a first-principle radiative transfer simulation 
and can be tested against observations. Here we examine a sample of 7 bursts with 
best GeV-|-optical data and test the model. We find that the observed light curves are 
in agreement with the theoretical predictions and allow us to measure three param¬ 
eters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, 
and the external density. With one possible exception of GRB 090510 (which is the 
only short burst in the sample) the ambient medium is consistent with a wind from 
a Wolf-Rayet progenitor. The wind density parameter A = pr'^ varies in the sample 
around 10^^ g cm“^. The initial Lorentz factor of the blast wave varies from 200 to 
540 and correlates with the burst luminosity. Radiative efficiency of the prompt emis¬ 
sion in the sample is between 0.1 and 0.8. For the two bursts with detected optical 
flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the ex¬ 
ternal blast wave. Remarkably, the model reproduces the entire optical light curve of 
GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV 
data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they 
can be detected by ground-based Gherenkov telescopes. 

Subject headings: plasmas - radiative transfer - plasmas - gamma-ray burst: general 


1. Introduction 


Cosmological gamma-ray bursts (GR Bs) radiate most of their energy in the soft g amma- 
ray band between 100 keV and 10 MeV ( Goldstein et al. 2012; Ackermann et al. 2013b). The 
MeV burst typically lasts seconds or minutes and is then followed by broad-band afterglow emis- 
sion, which is associate d with the deceleration of the explosion ejecta by the ambient medium 
(jMeszaros fc ReesI 119971 ). Afterglow observations can be used to estimate the main parameters of 
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the GRB explosion — its Lorentz factor, kinetic energy, and the density of the ambient medium. 
Interpretation of observations became, however, challenging after the Swift sat ellite discovered 
bizarre X-ray and optical light curves of the early afterglow (jGehrels et al.l |2009|) . Swift observa¬ 
tions challenge the standard assumption that afterglow is emitted by the decelerating shock wave 
in the external medium; instead, i t was proposed tha t afterglow is produced by a long-lived reverse 
shock (lUhm fc BeloborodovI 120071 : iGenet et al.l 120071 ) . Disentangling the two possible mechanisms 
is difficult, and a reliable method for deducing the explosion parameters from observations has been 
lacking. 

The discovery of GeV flashes by the Fermi satellite opens a new way for solving this problem. 
The observed flashes have similar light curves with a special shape: they begin with a delay and 
sharply peak well before t he end of the MeV burst; then the light curve exhibits a long gradual decay 
(Ackermann et al. 2013b). It is natural to associate the extended GeV emission with the external 


blast wave ( Zou et al.l l2009l : iKumar fc: Barniol DuranI l2009l : iGhisellini et al. 201C) although this 


scenario faced difficulties in explaining the early peak of the GeV flash (|Gao et al 


2011 


20001 : 


He et al 


Maxham et al.l[201ll ). The puzzling peak was recently explained by the exp onential loading 


of th e external medium by the prompt MeV radiation ahead of the blast wave (jBeloborodov et al 


20141 ) . Since the radius of pair loading is well dehned and can be determined from the prompt GRB 
observations, the GeV flash provides a standard “ruler” and a unique opportunity for disentangling 
the explosion parameters. Detailed m odeling of the pair-l oaded blast wave was performed for 


GRB 080916G ( Beloborodov et al.ll2014l ) and GRB 130427A ( Vurm et al.ll2014l ). It was shown that 
the GeV flash is emitted by the thermal plasma heated in the forward shock, which is cooled by 
inverse Gompton (IG) scattering of photons of lower energies. 

Along with the IG flash, the blast wave should emit synchrotron radiation, in particular in the 
optical band. Pair loading shapes the optical synchrotron light curve similarly to the IG GeV flash. 
As a result, the model predicts an optical flash with a sharp peak at a time close to the GeV peak. 
Following the peak, the optical emission should quickly decay, as the pair-loading effect is reduced, 
and the steep decay should be followed by a flatter light curve of the normal pair-free after glow. 
The predicted GeV-l-optical flash h as been detected i n GRB 130427A ( Vestrand et al. 2014 ). and 
the radiative transfer simulation of I Vurm et al.l (j2ni4l ) reproduced both the GeV and optical light 
curves of the flash. 


The consistency of the model with the data is remarkable given the fact that it has only four 
adjustable parameters: (1) the ambient wind density parameter A = pR^, (2) the explosion Lorentz 
factor Tej, (3) the prompt emission efficiency Crad = A^grb/(A^grb + Eej), which determines the 
isotropic kinetic energy of the explosion for a GRB with observed isotropic energy FIgrB) and 
(4) the magnetization of the blast wave cb- This fourth parameter only enters the calculation of 
the optical light curve and does not affect the GeV flash; therefore, the model of GeV emission has 
only three adjustable parameters. Other parameters of the explosion, e.g. the deceleration radius 
Rdec and the pair loading factor Z±{R), are not adjustable — they are calculated from the blast 
wave dynamics and the observed prompt radiation. 
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In this paper, we investigate all observed GRBs with well measured light curves of the flash. 
Our goal is to further test the proposed model and, if the model hts the data, to determine the 
parameters of the GRB explosions. We perform radiative transfer calculations individually for each 
burst, search for a blast wave model that would be able to reproduce the observed light curves 
of GeV emission and, if data is available, optical. Section 2 describes the sample of GRBs, and 
Section 3 describes the model and the method of our analysis. The results are presented in Section 4 
and further discussed in Section 5. 


2. GRB sample 


The choice of bursts for our sample is based on two criteria: (1) good Fermi LAT data, 
which provide the shape of the GeV flash, and (2) measured cosmological redshift z. In addition, 
we looked for bursts with available early optical data. We identified seven GRBs useful for our 
analysis: 080916G, 090510, 090902B, 090926A, 110731A, 120711A, and 130427A. 


The Fermi LAT data are taken from the published LAT catalogue (jAckermann et al.ll2013bl ) for 
GRBs 080916G, 090510, 090902B, 090926A, 110731A. Two recent GRBs 120711 A and 130427A are 


not in the catalogue; their publishe d LA T and optical flash dat a are taken from lAckermann et al 
(2014), Martin-Garrillo et al. (2014), and Vestrand et al. (2014). 


The prompt radiation data (which ar e used as an input of our transfer simu lations) are 


from 


Abdo et al.l ( 2009bl ) for GRB 080916C, lAckermann et al.l ( 2010l ) for GRB 090510, Abdo et al 


are taken from Golenetskii et al 


( 2009a ) for GRB 09 0902B, Ackermann et al. ( 201ll ) for GRB 090926A. Ackermann et alT ( 2013a ) for 
GRB 110731A, and Gruber &: Pelassal ( 2OI2I) for GRB 120711A. T he prompt data for GRB 130427A 


(j2013l ) ; lAckermann et al.l (j2014l ) . These papers provide an approx¬ 


imate description of the spectral evolution in each GRB, using spectral fits in temporal bins. This 
evolution is useful to take into account in our transfer simulations, however its details are not 
essential and weakly affect the results. For instance, in GRB 090510 it is safe to neglect the small 
precursor, as it carries small energy. In GRB 110731A we merged the first two time bins “A” 
an d “B” into a single bin to smoothen the jumps in luminosity and spectral parameters reported 
by lAckermann et al.l (j2013al ). These jumps mostly result from two different ways of fitting the 
observed spectrum: bin A spectrum was fitted by a power law with an exponential cutoff while bin 
B spectrum was fitted by a Band function whose high-energy slope is affected by the inclusion of 
the LAT data (see Ackermann et al. 2013a) 0 


The optical and X- ray afterglow data ar e from iGreiner et al.l (l2009l) for GRB 0809 16C, lDe Pasauale et al 


(201C) for GRB 090510, 

r*andev et al. 

(2010 for GRB 090902B, 

Swenson et al. 

(2010 for GRB 090926A, 

Ackermann et al. 

(2013? 

J) for GRB 110731A, 

Martin-Garrillo et al. 

(2014) for GRB 120711 A. and 


^The inclusion of GeV data in the Band component can be dangerous, as the GeV signal contains a separate 
component from the external blast wave, which can corrupt the inferred parameters of the Band MeV spectrum. 


































































































- 4 - 


Vestrand et alJ ( 20141 ) for GRB 130427A. We use the afterglow data to estimate the soft radia¬ 


tion field in the source, which dominates the IC cooling of the blast wave after the prompt MeV 
radiation. The afterglow reconstruction uses the simple power-law interpolation of the observed 
spectra and light curves. Extrapolation of the available data to earlier times was needed for several 
bursts. The disadvantage of this method is that the accuracy of the simple power-law extrapolation 
is uncertain, in particular for GRBs 090510, 080916G, 090902B, and 090926A. The advantage is 
that the method is well defined, model-independent, and minimizes special treatment of bursts 
with incomplete afterglow data. Fortunately, the rough reconstruction of the target soft radiation 
for IG scattering does not create large uncertainties in the light curve of the predicted GeV flash. 
In particular, in the fast cooling regime, the details of the target spectrum play almost no role 
( Beloborodov et al. 2014 ). 


3. The model 


Pair loading of the external medium by the MeV radiation can be accurately calculated for any 
observed GRB with a known redshift, using the observed luminosity (isotropic equivalent) and spec¬ 
trum of the prompt radiation. For any optically thin medium, the pair loadi ng factor Z+ = n+ln 
does r iot depend on the m edium density n, and is only a function of radius R (jThompson fc Madau 
20001 : iBeloborodovIl2002l ). The function Z±{R) is obtained by solving radiative transfer for the 


prompt MeV radiation. The transfer weakly affects the observed prompt radiation, however it 
strongly impacts the medium by depositing momentum and creating new particles. The pair load¬ 
ing factor is huge at small radii, Z± ~ 10^ — 10®, and is steeply reduced outside a characteristic 
“pair-loading” radius. The fast evolution of the density shapes t he peak of the GeV flash and 
its initial decay, as described in detail in [Beloborodov et al.l (|2014l ). For a typical bright burst 
observed by Fermi LAT, the peak radius Rp is comparable to 10^® cm. It is smaller than the 
deceleration radius of the blast wave, i?dec! and therefore the GeV flash peaks before the onset of 
normal afterglow, which is shaped by the blast wave deceleration. 


The mechanism of the GeV flash may be summarized as follows. The blast wave has a high 
Lorentz factor T and heats the pair-loaded external medium to a relativistic temperature. The 
medium is cooled behind the shock via inverse Gompton (IC) and synchrotron emission. IC cooling 
is extremely fast as long as the blast wave is exposed to the prompt MeV radiation (which is 
produced at much smaller radii and gradually overtakes the blast wave, with the relative speed of 
Au = c/2r^). When the prompt GRB radiation fully overtakes and decouples from the blast wave, 
the shock-heated plasma is cooled via the slower synchrotron-self-Compton (SSC) emission. The 
SSC regime occurs in the far tail of the GeV flash; the tail is well observed for some GRBs, in 
particular in GRB 130427A. 


In contrast to the older concept that G eV emission comes from nonthermal particles acceler¬ 
ated in the shock, iBeloborodov et al.l (120141 ) showed that the main, thermal population produces 
GeV-TeV emission. Therefore a large fraction of the blast wave energy is radiated in the high-energy 
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bands. The thermal plasma also produces the synchrotron optical flash. Electrons/positrons in¬ 
jected with the thermal Lorentz factor strongly dominate synchrotron spectrum at frequencies 
~ hy^j^ei^/meC, which covers the optical band during the flash. 

The dominance of radiation from the thermal plasma makes the model simple to test, without 
invoking phenomenological parameters describing the nonthermal tail. Nonthermal particles dom¬ 
inate synchrotron emission ai u > Um] this regime applies to the late optical afterglow when Vm 
decreases below 10^® Hz. During the flash, v > Um only at high (X-ray) frequencies. An extended 
nonthermal tail of the electron distribution can produce synchrotron radiation from the X-ray band 
up to ~ 0.1 GeV, contributing to the flux detected by LAT. This radiation is, however, found to be 
negligible in our best-fit models of GeV flashes (except perhaps the special case of GRB 090510) H 
Therefore, effectively three parameters (Tej, Cradj A) enter the fitting of the GeV flash, and is 
only relevant for the optical flash. 


Nonthermal synchrotron radiation implicitly enters our flash model in a different way: it 
provides targets for IC scattering in the SSC tail of the GeV flash. Since the nonthermal synchrotron 
modeling is expensive (and uncertain, especially when the reverse shock contribution is included), 
we estimate the targets using the actual observed afterglow (Section 2). 


The model must assume a val ue for the fraction of the shock energy that is given to the 
thermal electron/positron plasma. iBeloborodov et al.l (120141 ) showed that Ce ~ 1 during the peak 
of the flash. At later times (when Z± dro ps to 500), is reduced to 0.3, a typical value reported 


by the simulations of collisionless shocks (jSironi &: Spitkovskvl 120111 1. Note that the thermal is 
much less uncertain than the corresponding parameter of nonthermal particles; it is frozen and 
taken the same for all bursts. 


The progenitor wind has the mean molecular weight per electron /ig = p/neTUp = 2 (elements 
heavier than hyd rogen). The corre ct choice of pe is esse ntial in simulations of G RB afterglow (see 
the discussion in IVurm et al.ll2014l and comparison with iPanaitescu et al.ll2013l l. GRB 090510 is a 
short GRB and a special case; therefore for this burst we also search for a solution with uniform 
external medium and pe = 1- 

Appendix A summarizes simplified analytical estimates for the theoretical GeV-|-optical flash, 
which demonstrate basic trends in the model. Deviations of the accurate simulations from the 
estimates highlight the importance of detailed calculations for each GRB individually, using as an 
input its observed prompt radiation. Below the model with adjustable parameters A, Tej, Crad (and 
es, if optical flash data is available) is calculated for each GRB and fitted to t he data. The calcu¬ 


l ation involves a careful radiative transfer simulation, as described in detail in IBeloborodov et al 


(|2014l b The simulation is expensive and we do not attempt a formal fitting of the data that would 
give for the best fit. Instead, we manually search the parameter space for an acceptable solution. 


^During the peak of the flash, synchrotron cooling of the blast wave is negligible compared with IC cooling (as 
long as e_B ^ 0.1). It remains negligible in the tail if es <C 10“®. 
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Table 1: Three main observational parameters of the GRBs in our sample and four adjustable 
parameters of the model. T'grb ~ prompt radiation energy (isotropic equivalent), Tgrb = T 9 o/(l + 
z) - redshift-corrected duration of the prompt emission, z - cosmological redshift. A = pE? - wind 
density parameter, Tej - ejecta Lorentz factor, Crad ~ prompt efficiency, eb - magnetization. 


GRB 

7?GRB 
[10®"^ erg] 

Tgrb 

[s] 

z 

A 

[10^1 g/cm] 

Tej 

Crad 


080916C 

00 

00 

12 

4.35 

1.5 3.5 

900 1400 

0.17 

- 

090510 

0.11 

1.1 

0.903 

1.2 ^ 2 

700 800 

0.1 

- 

090510 (uniform) 



- 

n = 2 X 10^ cm“® 

900 

0.1 

- 

090902B 

3.6 

7.8 

1.822 

1 2 

600 -)• 900 

0.4 

- 

090926A 

2.2 

4.2 

2.106 

1 2 

600 1000 

0.25 

- 

110731A 

0.76 

1.9 

2.83 

0.4 0.8 

800 1100 

0.2 

- 

12071 lA 

1.65 

48 

1.405 

1 3 

320 -)• 400 

0.3 

10-5 ^ 10“® 

130427A 

0.85 

15 

0.34 

0.15 0.5 

300 350 

0.8 

10-5 ^ 2 X 10-"* 


Table 2: Other physical quantities calculated from the model: Tp ~ blast wave Lorentz factor 
at radius Rp of the GeV peak, R^ec ~ deceleration radius, R± - pair loading radius, tgc - time 
(measured in central engine frame) when the IG cooling of thermal electrons becomes inefficient. 


GRB 

Bp 

Rp 

[IQl® cm] 

-^dec 
[10^® cm] 

R± 

(lO^® cm] 

tsc 

[s] 

080916C 

540 

1 

10 ^ 6 

10 

200 -J- 300 

090510 

500 

0.2 

0.35 -f- 0.3 

0.9 

100 200 

090510 (uniform) 

400 

0.2 

0.3 

0.9 

> 10^ 

090902B 

300 

1 

4^3 

9 

6 X 10® 3 X 10'* 

090926A 

390 

1.1 

5^4 

9 

1.5 X 10® ->■ 6 X 10® 

110731A 

540 

0.5 

2.5 2 

4 

600 ->-4x 10® 

12071 lA 

200 

3 

10 <- 6 

12 

350 ^ 10® 

130427A 

250 

1 

3 2.5 

6 

2 X 10® -J- 7 X 10® 


4. Results 

For each GRB in the sample, the model successfully reproduced the observations in a nar¬ 
row region of the parameter space, which allowed us to measure the parameters (Tabled]). The 
theoretical light curves of the GeV and optical flashes are compared with observations in Figures 

For each GRB, we found two solutions whose predictions for the LAT/optical flux differ by 
about Icr. This gives a rough estimate of the “error bar” on the measured A and eb- A moderate 
uncertainty also results from the uncertain contribution of the prompt emission to the GeV light 
curve at early times, which is suggested by the variability in the observed light curve. It implies a 
somewhat lower peak of the smooth flash associated with the external blast wave. 

In all seven GRBs the GeV peak occurs while the prompt emission is still going on. This fact 
alone shows that the GeV flash is not associated with the deceleration radius of the blast wave. 
Table [2] gives more details of the reconstructed explosion — Lorentz factor Tp and radius Rp of the 
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blast wave at the GeV peak, the deceleration radius Rdec, the pair loading radius R± where Z± 
drops below 2, and the (redshift-corrected) time tgc when the cooling (synchrotron+IC) of thermal 
electrons becomes inefficient. We find Z±{Rp) ~ 10^ — 10^, a robust feature of our model: pair 
loading of the external medium is responsible fo r the extremely efficient p rocessing of shock energy 
into GeV gamma-rays, as explained in detail in iBeloborodov et al.l (|2014l i. 


The pair loading Z± rapidly decreases during the peak and early decay of the GeV light curve; 
as a result, the shock energy per lepton increases and the characteristic IC photon energy sweeps 
across the GeV band. All the presented models show similar evolution of the flash spectrum: the 
spectrum is soft during the rise of the GeV flash, then quickly hardens and remains approximately 
flat [vLy ~ constant) during the most luminous phase. This behavior is consistent with the 
spectral slopes observed by Fermi LAT (jAckermann et al.ll2ni3bl ) (with one notable exception — 
GRB 090510, whose spectrum is softer, see Section [53]) . At later times the predicted spectrum 
slightly hardens, approaching vLy oc that characterizes emission from fast-cooling thermal 
electrons. The predicted behavior of the optical spectrum is similar: flat near the peak, followed by 
moderate hardening as the synchrotron frequency Vm moves beyond the optical band. At present 
we possess no data on the spectral slope of the optical flash; its predicted behavior could be tested 
by future optical observations. 


Figures 1-3 also show the predicted TeV emission to provide some guidance for future obser¬ 
vations with Gherenkov telescopes. The characteristic IG photon energy reaches the TeV band 
in seconds or minutes after the GeV peak. The high-energy spectrum is most extended near the 
pair loading radius R±] the maximum photon energy is typically between a few hundred GeV and 
ten TeV, and scales as E\d-a/{A-E q^-q) where E\d-a is the kinetic energy of the blast wave. The 
electrons behind the shock are still fast-cooling at this stage, and their bolometric IG luminosity 
is ~ ee£'kin/(4t). The flat (or even slightly rising) high-energy spectrum places a sizable fraction 
of this luminosity into the TeV band, resulting in comparable GeV and TeV fluxes. Our transfer 
simulations also show that the high-energy flash is partially absorbed by local photon-photon col¬ 
lisions. This intrinsic absorption is included in all GeV-TeV flash models presented in this paper, 
however the effect is never dramatic, because the flash always peaks where the radiation column 
density is reduced below a critical value (see Beloborodov et al. 2014 for a detailed discussion). In 
particular, intrinsic absorption does not create a break of the high-energy spectrum, although it 
does moderately affect the spectrum shape. 

The TeV emission peaks immediately after its onset. Its decay is relatively slow, so that most 
of the energy above 0.1 TeV is emitted approximately over a decade in time after the peak, typically 
within a few minutes of the GRB trigger (Figures 1-3). The TeV flux from the thermal electrons 
behind the shock cuts off when the characteristic energy of the IC photons falls below 0.1 TeV, 
which typically occurs after a few hours or a day. An additional nonthermal electron population 
could affect the rise of the TeV light curve and extend the emission to later times, however it would 
not influence the most luminous phase. 








The seed photon field for IC scattering significantly changes at the time when the prompt 
MeV radiation completely overtakes the blast wave, as the remaining X-ray and optical afterglow 
radiation is softer and less luminous. The theoretical GeV light curves show no special features at 
this transition (Figures [T](3]) . This is because the shock-heated plasma is still in the fast-cooling 
regime; therefore it produces the same IC emission regardless of the details of the target spectrum. 


At the time of the GeV peak the ejecta kinetic energy is still being transferred to the blast 
wave via the reverse shock. Comparison of the ejecta and blast wave Lorentz factors in Tables [T] and 
[2] reveals that the reverse shock is at least mildly relativistic in most GRBs in the sample. When 
Tej ^ T (ultra-relativistic reverse shock), the blast wave dynamics and radiation are insensitive to 
Tej. Then the model of the GeV flash effectively has only two parameters: A and Cj-ad- The flash 
observations provide accurate estimates for the blast wave Lorentz factor T, however Tej is more 
difficult to measure. In some cases, only lower limits on Tej are obtained. 


The relativistic reverse shock crosses the main part of the ejecta (carrying most of the explosion 
energy) in about the same time as it takes the prompt MeV radiation to fully overtake the blast 
wave. This time also corresponds to the observed time Tgg and the deceleration radius i?dec- 
After this moment, the blast wave is approximately described by the Blandford-McKee self-similar 
solution. Note that i?dec < R± for all bursts in the sample, i.e. the effects of pair-loading continue to 
impact the afterglow emission for some time after the end of the prompt emission. This “memory” 
of pair creation is even more signi ficant at low frequencies where the emission occurs in the slow- 
cooling regime (jBeloborodovI 120051 ). 


Discussion 


In this paper, we tested the GeV flash model proposed by iBeloborodov et al.l (j2014l i using 
a sample of GRBs which includes all bursts with good LAT data and known redshifts. All 7 
bursts in the sample are intrinsically bright, with peak luminosities ranging from 1053 

erg s ^ 

(GRB 130427A) to ~ 10^^ erg s“^ (GRB 080916G). We performed radiative transfer simulations 
for each of the 7 bursts. The input of the simulation is the observed prompt MeV radiation and 
the observed optical/X-ray afterglow, and the output is the GeV light curve. The model has three 
adjustable parameters: the Lorentz factor of the GRB ejecta Tej, the ambient density parameter A, 
and the radiative efficiency of the prompt emission Crad = -Egrb/(T'GRB + -®ej)- We found that the 
model well explains the observed light curves in the sample. This allowed us to obtain estimates 
for Tej, A, and e^ad- 


5.1. Ambient medium 

Explosion into a wind-type medium is consistent with the observed flash for all GRBs in the 
sample except possibly GRB 090510 (see below). We found that the wind parameter A shows 
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moderate variations in the sample, between 0.15 x 10^^ g cm“^ and 3.5 x 10^^ g cm“^ (Table 1). 
These values are compara ble with typical ^ ~ 3 x 10^^ g cm“^ estimated for the winds of Wolf- 
Rayet stars in our Galaxy (|Crowtheiil2007l i . This result provides further support for the association 
of GRBs with collapse of Wolf-Rayet stars. Evidence for this associatio n was previously provid ed 
by a few GRBs with a detected supernova counterpart of type Ib or Ic (jWooslev &: BloomI 120061 1. 


We also note th at the obtained values of A do not contradict the upper limits estimated by 
Hascoet et al.l (|2014l i with a different method and for a different sample of bursts. In that work, 
the most constraining upper limits ^max ~ 10^^ g cm“^ were obtained for GRBs with luminosities 
below 10^^ erg s“^ (Figure |3|). Comparison of these upper limits with our estimates presented here 
suggests that the luminous bursts detected by LAT have systematically higher A. 


In most bursts in the sample the flash peaks much earlier than the blast wave reaches the 
deceleration radius Rdec- In this situation, the models of uniform and wind external medium 
predict very different light curves of the flash, and we find that only the wind medium is consistent 
with the data. The distinction is less clear when Rp is comparable to Rdec- Then both the shape 
of the GeV peak and the decay of the light curve are relatively insensitive to the profile of the 
ambient medium (provided that the density at Rp remains similar and the emitting electrons are 
fast-cooling). In our sample, GRB 090510 falls into this category (see below). The separation 
between Rdec and Rp is also moderate in GRB 120711A and GRB 130427A, however the wind 
medium is still preferred in these two cases. Thus the wind medium is preferred by the analysis of 
all long bursts in the sample. 


5.2. Lorentz factor 

Measurement of GRB Lorentz factors is a long standing problem. Until recent work, the main 
method was based on estimates of photon-photon opacity for the gamma-rays detected by LAT , 


which gave lower limits on Tej (e.g. iLithwick &: Sari 


2001 


Granot et al 


2008 


Hascoet et al 


201J). 


Reconstruction of the GeV hash mechanism gives a valuable measurement of the Lorentz 
factor of the blast wave T before its deceleration, and also provides an esti mate for the ejecta 


2014 

ej ~ . .... 

) and GRB 130427A ( 

Vurm et al. 

2014) 


7 GRBs. Figure [5] shows the obtained Lorentz factors versus the burst luminosity. We observe 
a positive correlation between T and the average luminosity of the GRB Rgrb = Rgrb/Tgrb, 
where the prompt burst energy Rgrb and its approximate duration Tgrb are given in Table 2. 
Provided that Tej is not much higher than T (which is the case for GRB 130427A, and likely in the 
other bursts) we find a LGRB-Bej correlation, which may be roughly approximated by a power-law 
relation Tej ~ . Future observations of GeV hashes may allow a better measurement of this 

correlation. 


There is another method of estimating T using the peak time of the optical afterglow (ILiang et al. 




























- 10 - 


2010 : Ghirlanda et alJl2012 ). The method is based on the assumption that the peak is emitted at the 


deceleration radius of the blast wave. This assumption can be invalid for bright bursts which have 
large pair-loading radii (as demonstrated for the sample of GRBs studied in this pap er); however it 


may b e reasonable for less luminous bursts. A recent refinement of this method bv iHascoet et al 


(|2014l ) gave measurements and upper limits for T in a large sample of bursts, which are included in 


Figure Ho They extend the Lgrb-B diagram to lower luminosities Tgrb < 10^^ erg s“^, however 
do not allow a reliable measurement of the correlation in this region because of the large number 
of lower limits. Overall, the data is consistent with the existence of a lower bound T > Tmin which 
increases with Tgrb- This bo und may result from strong subphotospheric adiabatic cooling in 
explosions with Tej < Tmin Isee lHascoet et al.ll2014l ). 


We also note that within uncertainties our estimate for the jet Lorentz factor in GRB 130427A 
is in agreement with the value Tej = 450 found from radiative transfer modeling of its prompt 
emission (Vurm & Beloborodov, in preparation). 


5.3. Radiative efficiency 


The inferred radiative efficiency Crad of the prompt GRB emission varies between 0.1 and 0.8 
in the sample, i.e. Fiej/R^GRB = 0.25 — 9. Such high values of radiativ e efficiency may be expected. 
High Crad was previously suggested by the late afterglow a nalysis le.g. iRacusin et al.ll201ll) and also 


expected in theoretical models of the prompt emission (cf. BeloborodovIboiC ; Vurm & Beloborodov, 
in preparation). 


5.4. The special case of GRB 090510 

The fG cooling of shock-heated wind medium well explains the observed GeV flash for all 
GRBs in our sample except one: GRB 090510. In this case, the model rather well reproduces the 
peak of the flash at tobs ^ 1 s and its initial decay at Gbs < 10 s, however at Gbs > 20 s the 
theoretical emission falls short of the observed flux. 

GRB 090510 is also special as the only short burst in our sample. Short GRBs are normally 
not associated with massive progenitors, although some of them may be “impostors” in the short 
class (the impostors are bursts that have a massive progenitor but happen to have short duration). 

The dehciency of the theoretical GeV emission in GRB 090510 may be explained in two ways: 

(1) The ambient medium is not a wind from a massive progenitor. Indeed, we hnd that the 


® This method gives estimates for where pdec is the density of the external medium at the deceleration 

radius. Therefore, there is a weak dependence of the inferred F on the assumed pdec- In the left panel of Figure 4 we 
used Pdec ~ g cm“^ which corresponds to T 10^^ g cm“^ and Rdec ~ 10^® cm. 
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entire light curve of the GeV flash is reasonably well reproduced if the ambient density is uniform, 
p ~ const, rather than wind-like, p oc The uniform density must, however, be quite high, 

n = p/nip ~ 2 X 10^ cm~^, well above the typical density of interstellar medium. 

(2) An additional emission component is present in the flash observed by LAT. It can be 
synchrotron emission (extending to ill > 100 MeV) from nonthermal electrons accelerated in the 
blast wave. This interpretation is consistent with the relatively soft photon index (/3 ~ 2.5 — 3) 
measured for GRB 090510 at the time when the additional component dominates. Our simulation 
did not include possible nonthermal electrons because they are harder to model form hrst principles 
and require additional phenomenological parameters. We also note that the nonthermal component 
is not needed for the other 6 bursts — their GeV flashes are well explained by pure IG emission 
from the thermal plasma. 


5.5. Optical flash 


Two GRBs have a well measured optical counterpart of the GeV flash. The expected light 
curve of this (synchrotron) counterpart is obtained directly from the model of the GeV flash by 
introducing one additional parameter eg- 


GRB 130427A has the best coverage in the LAT and optical bands. The entire GeV flash 
(f < 1 d), and the main peak and steep decay of the optical flash [t < 100 s), are well reproduced 
by the model of the pair-loaded forward shock. The optical light curve requires a n additional 
contr ibution after 100 s, which we interpreted as emission from the reverse shock (jVurm et al 


20141 ). The forward shock produces the light curve similar to that observed in GRB 120711A 


(discussed below), with a plateau ending around 10^ s. This plateau may be responsible for the 
hump in the optical light curve of GRB 130427A at f ~ 10^ s. 


The result obtained for GRB 120711A is rather striking. The model reproduces the entire 
complicated optical light curve: the sharp rise at 30 s, the peak at 50 s, the steep decay between 50 
and 300 s, the plateau and the break at lO'^ s. Our model has only four adjustable parameters A, 
Crad; Tej, and cb- Such a remarkable ht of the optical light curve can hardly be obtained by chance. 
The same model also hts well the decay of the GeV flash, both the slope and the normalization. 
The GeV data are less constraining in this burst as the peak of the GeV flash was missed by LAT 
observations. 


The inferred ~ 4 x 10“^ for GRB 130427A and eg ~ 3 x 10“® for GRB 120711A give rough 
but reliable estimates for the characteristic magnetization in the externa l blas t wave. In a more 
detailed model cb may be changing with time as discussed in IVurm et al.l ()2014l ). The excellent ht 
for GRB 120711A does not require such variation. 







- 12 - 


5.6. TeV emission 


Our model predicts luminous TeV emission accompanying the GeV flash, which can be detected 
by ground-based Cherenkov telescopes. The bulk of the TeV fluence is accumulated within 1-10 
minutes after the GRB trigger; the energy radiated above 0.1 TeV in the sample ranges from 10^^ 
to 4 X 10^^ erg, and constitutes up to 30% of the prompt MeV fluence. 

The predicted efficiency of TeV emission varies substantially from burst to burst. For exam¬ 
ple, in GRB 120711A most of the energy above 100 MeV is radiated in the TeV band (Figure 
[3|). In contrast, TeV emission is weak in the high-A (dense wind) models for GRB 090902B and 
GRB 090926A (Figured]). The difference in the TeV flux arises mainly from the different maximal 
IC photon energy that the thermal electrons can produce, Fimax ~ TythnieC^, where 7th is the 
Lorentz factor of the thermal electrons behind the shock. In GRBs exploding into denser winds the 
blast wave is slower and loses a larger fraction of its kinetic energy at an early radiative stage. The 
reduced F gives a lower Fimax so that it may not exceed a few 100 GeV (in the GRB rest frame). 

Nonthermal electrons can extend the IC emission above Fimax; this emission was not included in 
our model, and thus the 0.1 TeV light curves shown in the figures should be viewed as a lower limit. 
The shock energy given to nonthermal electrons is small compared with the thermal population. 
Therefore, their presence in the 0.1 TeV light curve becomes important only when the contribution 
from the thermal electrons is suppressed, i.e. when Fimax <0.1 TeV. 


A major factor limiting the GRB detectability in the TeV band is the extinction by extra- 
galactic background light (EBL). F or example, at redshift z = 1 the attenuatio n factor is ~ 0.5 at 


0.1 TeV and ~ 5 x 10“^ at 0.3 TeV (jPommguez et al.ll2011l : lGilmore et al.ll2009lf . This leaves a nar¬ 
row window between the Cherenkov detector threshold (presently ~ 50 — 100 GeV) and ~ 0.2 — 0.3 
TeV for most bursts, except those that happen at unusually small 2 ; (such as GRB 130427A). As 
long as the intrinsic high-energy spectral cutoff is well above 100 GeV, the observed turnover would 
arise from the extragalactic absorption and could be used to place independent constraints on the 
EBL density. 

The (approximate) low-energy threshold ~ 50 GeV, repositioning time of a few tens of seconds 
and sensiti vity of the currently operating Imaging Atmospheric Cherenkov telescopes (lACT) such 


as MAGIC (Aleksic et al 


2012; 


Sitarek et al. 


H.E.S.S. ( Hinton &: the HESS Collaboration 


2013), VERITAS (iHolder 


2011 


Kieda et al.ll2013l i and 


2004) 1 would have been sufficient to detect at least two 


bursts in our sample: GRB 120711A and GRB 130427A. Their predicted fluxes above 0.1 TeV a 
few minutes after the GRB trigger are well above the sensitivity limit despite the substantial EBL 
absorption for GRB 120711A {z = 1.405). In the case of GRB 130427A {z = 0.34) the TeV emission 
remains at a detectable level for its entire duration, i.e. until the cutoff at a fewx 10^ s. The cuto ff 


is consistent with the upper limit provided by the VERITAS observation at 1 day (lAliu et al.ll2014l ) 


The projected sensi tivity of the next generation Cherenkov Telescope Array (CTA) (|Eunk et al 


2 OI 3 I : llnoue et al.l 120131 ) is sufficient to detect most of the bursts in our sample (including GRB 
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090510 if it exploded into a wind medium), with the possible exception of GRB 080916C and 
GRB 110731A owing to their high redshifts. Gurrent efforts to reduce the energy threshold of the 
Gherenkov telescopes are key for future routine detection of high-energy GRB emission from the 
ground. 


In this paper, we focused on the flash observations, and only in the bands where we think 
it is dominated by the thermal plasma behind the forward shock. The analysis of all afterglow 
observations, at all times and energies, from radio (Laskar et al. 2013) to hard X-rays (Kouveliotou 
et al. 2013) is deferred to future work. It will be significantly more involved, as it has to include 
the emission from nonthermal particles, from b oth reverse and for ward shocks. The prospects for 
such a model for GRB 130427A are outlined in IVurm et al.l (120141 1: they argue that the proposed 
blast wave model can be consistent with radio and hard X-ray d ata with reasonable assumptions 
regarding the reverse shock and nonthermal particle acceleration. IVurm et al.l (120141 ) also pointed 
out the importance of using the correct mean molecular weight per electron, //g = 2 , expected for 
a Wolf-Rayet wind. The external density estimated from the late nonthermal optical flux scales as 
fj-l- Using /ie = 2 instead of /ig = 1 increases the inferred parameter A by the factor of 8 and, at 
least in the case of GRB 130427A, makes it consistent with the value measured from the GeV flash. 


We are grateful to Nicola Omodei for providing the LAT catalogue data, Antonio Martin- 
Carrillo for providing GRB 120711A data, and Tom Vestrand for providing GRB 130427A data. 
This work was supported by NSF grant AST-1412485, NASA ATP grant NNX15AE26G, and Swift 
Cycle 10 grant NNX14AI94G. 


A. Analytical estimates 


Here we summarize the analytical estimates derived in Beloborodov et al. (2014) to show basic 
trends in the flash model. Consider an idealized model of the prompt MeV radiation with a fixed 
spectrum and duration. One main parameter is left free — the prompt luminosity Tgrb- 


The Lorentz factor and radius of the bl ast wave at the GeV pea k are related to the observed 
peak time Tp by Equations (48) and (49) in iBeloborodov et al.l (120141 1. 


/ 'T \ ~ 3/13 

G «= 500 lIP (^) . (Al) 

/ rp Ijn 

Rp ^ 10^6 l 6/13 ^ j (A2) 

where L 54 = Lgrb/ 10^^ erg s“^ and Tp is measured in seconds. Note that the external density 
parameter A does not enter these relations; the values of Rp and Tp are controlled by pair loading 
Z±{R) which does not depend on density. The results of detailed simulations reported in Section 4 
for our GRB sample show moderate deviations from the simplified relations (|Aip and (|A2p . 
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The photon number emitted in the GeV flash (its isotropic equivalent) is estimated as 


NceV ~ 


4:7rARp 

flelTlp 


^±{Rp) -M. , 


(A3) 


where 


M{E) ~ 




(A4) 


is the average number of IC photons of energy E' ~ 1 GeV emitted by a single post-shock electron, 
and Et ^ 1 MeV is the typical energy of target prompt photons. Combining Equations (IAip - (IA4l) . 
one can express the wind density parameter as 


All 


0.3 iVceV,56^54*'''" 


\ 1 

-4/13 

Z±{Rp) 

l + z_ 


10^ 


Et y/V E \ 

1 MeV J Vl GeVy ’ 


(A5) 


where we have normalized the number of high-energy photons to the typical value Nq^v r\j 1056 
observed in our sample. 


The kinetic luminosity of the ejecta Lgj can be roughly estimated ass uming a relativistic revers e 
shock and a pressure balance between the reverse and forward shocks fsee lBeloborodov et al.ll2014l l. 


Lgj ~ 4 X 10^^ A'GeV,56 


l + z 


-1 


Et \ 
1 MeV ) 


1/2 


E \ 


1/2 


1 GeV 




erg s 


-1 


(A6) 


If the reverse shock is not ultra-relativistic, equation 


should be considered as a lower limit. 


Finally, if the simultaneous opt ical flash is detected, an estimate of the blast wave magnetiza¬ 
tion is given by (see Equation 81 in iBeloborodov et al.ll2ni4l ) 


eu ~ 3 X 10 ^ 


E±{Ropt) 

100 


L 


2 

opt,49 


R 


-2 

Opt, 16 


^11 



Et 

1 MeV 


(1 + ^) 


(A7) 


where Topt ,49 is the peak luminosity of the optical flash normalized to 10^® erg s and i?opt is the 
radius where the optical flash peaks, which is slightly outside the peak radius of the GeV flash, Rp. 
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Fig. 1.— Model light curves and data for GRBs 080916C, 090902B, 090926A and 110731A. 
Theoretical photon flux (isotropic equivalent) and data above 100 MeV are shown in blue, and the 
photon flux predicted above 100 GeV is shown in magenta. Time t is measured in the rest frame 
of the central engine, t = tobs/(l + z). Solid cnrves show the high-A (high-Tej) model and dashed 
cnrves show the low-A (low-Tej) model (see Table 1). 
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Fig. 2.— Model light curves and data for GRB 090510. Theoretical photon flux (isotropic equiv¬ 
alent) and data above 100 MeV are shown in blue, and the photon flux predicted above 100 GeV 
is shown in magenta. Time t is measured in the rest frame of the central engine, t = tobs/(l + ^)- 
Solid curves show the high-A (high-Tej) model and dashed curves show the low-A (low-Tej) model 
(see Table 1). The dashed-dotted curve shows the model with uniform external medium. 
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Fig. 3.— Model light curves and data for GRBs 120711A and 130427A. Solid curves show the 
high-A (high-Fej) model and dashed curves show the low-A (low-Fej) model (see Table 1). Upper 
panels: theoretical luminosity and data above 100 MeV (blue) and the luminosity predicted above 
100 GeV (magenta). Lower panels: theoretical luminosity and data at 2 eV. 
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Fig. 4.— Left panel; blast wave Lorentz factor (before deceleration) versus the prompt GRB 
luminosity Lgrb = F'grb/Fgrb- The 7 LAX bursts studied in this paper are plotted as magenta 
filled squares. We included the short GRB 090510 and modeled its GeV flash with both wind 
(magenta square) and uniform (blue triangle) external medium. For comparison, we also show the 
estimates (open circl es) and lower limits (arrows) for F obtained for a sample of weaker bursts with 
a different method (jHascoet et al.M2014l ). Right panel: distribution of wind parameters in units 
of 10^^ g cm“^. The LAX burst s stud ied in this paper are plotted in magenta. Arrows show the 
upper-limits from iHascoet et al.l (120141 ). 












